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ABSTRACT 

Jr 

A  linear  theory  for  a  gyrotroryfBWO;  has  been  developed  and  is  presented 
in  this  report.  The  theory  solves  a  reduced  one-dimensional  Maxwell-Vlasov 
equation  in  the  form  of  a  linear  integro-differential  equation  using  Laplace 
-transformation.  The  relative  amplitudes  among  the  waveguide  modes  and  beam 
modes  are  completely  determined  and  enable  us  to  determine  BWO  oscillation 
conditions . 

A  design  is  made  based  upon  this  analysis  including  velocity  spread 
effects.  A  several  tens  kilowatt  power  level  gyrotron  BWO  may  be  possible 
and  find  useful  applications  as  rf  heating  sources  in  fusion  devices  and 
as  driving  sources  in  high  power  gyrotron  amplifiers^  , -7  rn~For — 
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1.  INTRODUCTION 

The  gyrotron  has  demonstrated  unprecedented  capability  of  producing 

high  power  microwaves  with  high  efficiency  at  high  frequencies.^  Most  gyro- 

trons  use  a  cyclotron  maser  instability  (CMI)  in  a  forward  wave  region 

requiring  external  feedback  and  thus  cavity  structure.  The  rf  field  in  the 

cavity  is  well  approximated  by  the  cavity  structure  alone  and  thus  the 

2 

theory  is  fairly  well  understood. 

A  gyrotron  forward  wave  amplifier  operates  in  a  waveguide  without  exter¬ 
nal  feedback.  However,  its  gain  (growth  rate)  is  essentially  determined 
by  one  growing  mode  which  can  be  determined  by  a  complex  dispersion  rela¬ 
tion  alone,  besides  some  uncertainty  in  the  initial  coupling  of  the  beam 
mode-to-signal.  Up  to  this  uncertainty  the  amplifier  theory  is  also  fairly 

3 

well  understood. 

A  gyrotron  backward  wave  oscillator  (BWO)  operates  also  in  a  waveguide 
but  the  electron  beam  itself  plays  the  role  of  feedback.  Due  to  its  self- 
feedback  mechanism,  a  gyrotron  BWO  does  not  require  an  external  feedback 
structure  and  thus  allows  us  to  tune  the  frequency  electronically.  In  this 
case,  at  least  all  three  backward  propagating  modes  (two  beam  modes  and  one 
waveguide  mode)  play  equally  important  roles  and  thus  their  amplitude  infor- 

4 

mat ion  is  essential  as  well  as  the  complex  dispersion  relation.  Previously, 
we  have  developed  a  linear  theory  which  enables  us  to  calculate  the  ampli¬ 
tudes  from  boundary  conditions.  This  solution  of  a  boundary  value  problem 
using  Laplace  transformation  allows  us  to  investigate  the  backward  wave 
oscillator  which  we  will  present  in  this  report. 

Conventional  BWO  requires  a  slow  wave  structure  (usually,  helical  wire) 
which  severely  limits  the  power  handling  capability.  Typically,  a  sweep 
oscillator  which  uses  a  conventional  BWO  delivers  a  power  level  of  a  few 


tens  milliwatt  at  35  GHz 


A  gyrotron  BWO  dots  not  require  any  alow  wave  structure  and  thus  its 
power  handling  capability  is  no  worse  than  the  gyrotron  Itself  (sub-mega- 
watt  level).  Actually  a  gyrotron  BWO  operates  far  above  cutoff  while  a 
gyrotron  operates  near  the  cutoff  and  therefore  the  wall  loss  of  the  gyro¬ 
tron  BWO  should  be  even  smaller  than  the  gyrotron.  Also  it  is  tunable 
both  by  voltage  and  magnetic  field  while  conventional  BWO  is  tunable  only 
by  voltage. 

A  backward  wave  oscillator  is  expected  to  be  intrinsically  less 
efficient  than  a  forward  wave  oscillator.  This  is  due  to  the  less  favor¬ 
able  operating  conditions  with  beam  bunching  in  a  strong  field  region  and 
power  extracting  in  a  weak  field  region.  This  handicap  is  compensated  by 
the  wide  band  electrical  tunability.  Therefore,  emphasis  should  be  placed 
on  broad  band  tuning  capability  with  wide  band  coupler  design.  This  study 
shows  that  a  few  tens  killowatt  level  of  Ku-band  gyrotron  tunable  over 
almost  an  octave  range  may  be  possible  with  an  existing  electron  gun.  Such 
a  device  may  find  useful  applications  as  a  driving  source  to  a  gyrotron 
amplifier  or  as  a  source  in  the  study  of  rf  heating  of  fusion  devices. 

As  a  source  of  gyrotron  forward  wave  amplifiers ,  one  may  recall  that 
a  high  power  amplifier  usually  operates  at  a  low  gain  regime  requiring 
moderately  high  power  source  to  drive  it.  As  a  source  of  rf  heating  of 
fusion  devices  one  may  recall  that  heating  characteristics  are  highly 
dependent  on  the  magnetic  field,  temperature,  density  of  the  fusion  device 
as  well  as  the  polarization  of  the  rf  field.  The  tuning  capability  of 
gyrotron  BWO  may  be  invaluable  to  tune  the  frequency  according  to  the 
operating  condition  of  the  fusion  device. 
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In  this  report  we  will  present  a  linear  theory  of  a  gyrotron  backward 
wave  oscillator.  The  beam  is  assumed  to  be  a  hollow  gyrating  beam  from  a 
typical  MIG  gun ,  and  the  waveguide  is  assumed  to  be  a  rectangular  waveguide. 
Actually,  a  cylindrical  waveguide  has  also  been  studied  and  the  formulation 
can  be  carried  out  in  such  a  way  that  one-to-one  comparison  can  be  made  all 
along  the  calculation.  The  rectangular  TE**  mode  is  chosen  because  of  its 
ease  of  mode  control  (no  mode  competetion  up  to  an  octave)  and  design  of  a 
wide  band  coupler.  A  test  model  of  a  coupler  has  been  made  and  demonstra¬ 
ted  wide  useful  band  coupling  capability  within  a  few  dB  coupling  loss  over 
an  octave  besides  a  couple  of  sharp  lossy  spikes. 

In  Section  2,  we  reduce  Maxwell  equations  with  an  electron  beam  as  a 
source  Into  one  dimensional  equations  governing  the  axial  growth  behavior. 
The  reduction  is  possible  under  an  assumption  that  a  single  waveguide  mode 
Interacts  resonantly  with  the  beam  and  the  result  of  Interaction  is  mainly 
growing  or  decaying  of  the  mode  axially  with  little  modification  in  the 
-  transverse  direction.  Laplace  transformation  is  performed  to  convert  the 
differential  equations  into  algebraic  equations  with  the  boundary  values 
Included  explicitly. 

In  Section  3,  we  calculate  the  perturbed  electron  beam  function  from 
the  linearized  Vlasov  equations.  With  harmonic  expansion,  the  beam  modula¬ 
tion  can  be  written  as  hysteresis  integral  of  the  form  (3.9)  linearly  pro¬ 
portional  to  the  rf  field.  Laplace  transformation  again  helps  to  reduce 
this  Integral  into  a  product  of  two  Laplace  transformed  functions. 

In  Section  4,  the  phase  space  integral  is  carried  out  to  obtain  the 
induced  current  and  thus  the  source  term  in  the  Maxwell  equation.  The  inte¬ 
gration  is  arranged  in  such  a  way  that  all  the  phase  space  integral  can  be 
done  for  a  "cold"  beam  and  the  velocity  spread  effect  can  be  introduced 


1.4 


aftczvard  through  "statistical  average"  procedure. 

In  Section  5,  the  source  term  is  combined  with  Maxwell  equations  and 
che  resultant  algebraic  equations  are  solved  and  the  Inverse  Laplace  trans¬ 
formation  is  performed  to  obtain  the  field  as  a  function  of  z.  The  complex 
dispersion  relation  is  identified  as  a  denominator  of  the .integrand  of  the 
inverse  Laplace  transformation  whose  zeros  determine  poles.  The  residues 
at  the  poles  determine  the  amplitudes  of  each  mode.  For  a  cold  beam,  the 
denominator  is  a  quartic  function  giving  four  poles  -  two  waveguide  modes 
and  two  beam  modes. 

In  Section  6,  the  dispersion  relation  is  analyzed  in  terms  of  scale 
variables  which  enable  us  to  compare  the  relative  strength  of  beam  coupling. 
Eliminating  the  weakly  coupled  oppositely  propagating  mode  ,  one  arrives 
at  a  simple  approximate  cubic  dispersion  relation  which  can  be  easily 
analyzed  analytically.  This  formulation  enables  us  to  compare  directly 
with  the  conventional  Pierce  type  approach  in  a  conventional  microwave 
device. . 

In  Section  7,  the  velocity  spread  effect  is  examined  with  a  "generalized 
Lorenzian  distribution"  function  (7.1).  Other  than  the  regular  Lorenzian 
distribution  function,  the  number  of  modes  are  more  than  four.  A  more 
realistic  distribution  function  with  n»2 ,  the  total  number  of  modes  are  six 
and  this  distribution  function  was  used  to  investigate  the  velocity  spread 
effect  in  the  design. 
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II.  MAXWELL  EQUATIONS  WITH  A  SOURCE  TERM 

In  this  section,  we  will  derive  one-dimensional  equations  from  Max¬ 
well  equations  governing  RF-field  growing  or  decaying  in  presence  of  a 
source  (an  electron  beam)  interacting  predominantly  with  a  single  mode, 
particularly  TE^  modes  in  a  rectangular  waveguit  Of  course,  the  induced 
current  of  the  beam  -  which  will  be  calculated  frou.  Vlasov  equation  in  the 
following  two  sections  -  is  modulated  due  to  the  coupling  with  RF-fleld 
and,  thus,  Maxwell-Vlasov  equation  forms  a  self-consistent  integro-differ- 
ential  equation  which  will  be  solved  in  Section  V.  This  type  of  single 
mode  analysis  Is  valid  on  the  basis  of  the  following  assumptions. 

In  general,  the  presence  of  the  beam  in  a  waveguide  alter  the  field 
pattern  due  to  its  space  charge  and  make  them  grow  or  decay  due  to  its 
Interaction.  Expanded  in  terms  of  empty  waveguide  modes,  since  they  form 
a  complete  orthonormal  set,  the  space  charge  effect  appears  as  a  mode- 
coupling  and  the  interaction  makes  them  grow  or  decay.  In  some  cases, 

-this  mode-coupling  is  important  in  studying  the  mode  competetion  which  will 
be  dealt  with  somewhere  else. 

However,  if  the  beam  is  tenuous,  only  those  modes  interacting  reson¬ 
antly  with  the  beam  will  participate  significantly.  In  practice,  wave¬ 
guide  modes  are  fairly  well-separated  and  the  beam  parameters  can  be  chosen 
to  be  resonant  only  with  one  Interesting  mode.  In  this  case,  the  effect 
of  mode  coupling  is  negligible  and  the  dominant  effect  will  be  the  growth 
or  decay  of  the  resonant  mode  as  s  result  of  interaction,  allowing  us  a 
single  mode  analysis. 

The  above  observations  enable  us  to  write  an  ansatz  for  the  stationary 

P 

fields  of  TEno  mode  in  a  rectangular  waveguide  shown  in  Fig.  2.1  as 


Ey(x,t)  -  -  e  Ey(z)ey(x) 

Hx(x,t)  -  e"iutHx(z)hx(x) 

H  (x,t)  »  -  ie”1WtH  (z)h  (x) 

Z  Z  Z 

(E  -  E  -  H  »  0) 
x  z  y 

where  ey(x),  hx(x),  and  h^(x)  are  TE®q  empty  waveguide  mode  given  by 

h  (x)  ■  cos  k.x 
z  **• 


„  ,  v  m/c  d  .  .  N  tu/c 
e  (x)  =  —hr~  —  h  (x)  -  - 
y  k2  dx  z  ka 


sink.x 


e  (x)  -  e  (x)  ■  -  £“■  sink  x 
x  w/c  y  kA 

2  o)^  2 

<ki  -  ~  k„) 

c 

and  k„  is  a  wave  number  of  TE®q  mode  in  the  empty  wave  guide.  Note  that 
the  z-dependencies  of  the  fields  in  Eq.  (2.1)  are  left  undetermined  func 
tions  to  allow  growth  or  decay  which  can  be  determined  self-consistently 
-with  coupling  to  the  source  term. 

i  a 

Substituting  the  ansatz  (2.1)  into  Maxwell  equations,  V  x  E  =  - — 

1  |  ^  c  3 
and  V  *  H  »  -  —  +  —  1,  and  using  (2.2),  one  obtains  a  set  of  one  dimen 
cot  c 


sional  equations. 


Ey(z)  -  ik„Hx(z) 


Ey(z)  -  Hz(z) 


Ey(z)  +  k*  Ey(z)J  ey(x)  -  eiutJy(*,t) 


Using  TE  mode  projection  operator,  one  can  write  the  last  equation  as 


2.3 


¥  (z)  +  k^¥  (z)  -  i  —  -  elutT  (z)/C 
dz2  y  y  cc  y  n 


(2.4) 


where 


V*>  E/0dXjody  •y<*>Jy<*.t> 

/*a  /"b  2  2 

c»  i/odxV5’  ^  r 


(2.5) 


Equations  (2.3)  and  (2.4)  form  a  Maxwell  equation  part  of  Maxwell-Vlasov 
equation,  governing  the  growth  or  decay  of  KF  field  with  a  source  term 
given  in  the  form  of  hysteresis  integral  which  will  be  calculated  in  the 
following  two  sections. 

This  type  of  one-dimensional  integro-differential  equation  is  best 
solved  hy  Laplace  transformation  defined  by 


t (k)  £  J q  dz  e“ikz  F(z) 
(Ink  «  0) 


(2.6) 


where  k  is  a  complex  variable  with  sufficiently  large  negative  imaginary 
value  to  guarantee  the  integral  in  (2.6)  to  be  well-defined  even  when  F(z) 
is  an  exponentially  growing  function  as  in  the  case  of  instability  due  to 
a  resonant  interaction.  The  Laplace  transformation  of  Eq.  (2.3)  and  (2.4) 


k„H  (k)  -  kE  (k)  +  iE  (0) 
*  y  y 

Hz(k)  -  Ey(k) 


(2.7) 


where 


(k2  -  k2)Ey(k)  -  P(k)  -  ikEy(0)  -  |f*  z-() 


P(k)  -  -  i  ~  ?  eiUt  Jy(k)/Cn 


(2.8) 


(2.9) 


and  J  (k)  is  the  Laplace  transformation  of  Jy(s)  which  is  defined  by  (2.5). 

Mote  that  the  boundary  conditions  Ey(0)  and  (dEy/d z>2-0  "  co“* 

in  naturally  in  the  above  equations  which  might  have  been  lost  if  we  had 
used  Fourier  transformation  instead  of  Laplace  transformation.  Furthermore 
Fourier  transformation  similar  to  (2.6),  while  the  integration  runs  from 
—  to  +«,  cannot  be  made  well-defined  in  case  of  instability.  In  the  follow- 

will  calculate  the  source  term  P(k)  defined  by  (2.9) 


ing  two  sections,  we 
from  Vlasov  equation. 
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III.  VLASOV  EQUATION  -  PERTURBED  ELECTRON  DISTRIBUTION  FUNCTION 
In  a  linear  approximation,  the  perturbed  part  of  electron  distribution 
function  la  given  by 

-*•  ■? »  -+  5  f 


r  t 

^(x.p.t)  -  e  I 

)  t-s 


dt*  (E*  +  *■ —  x  H') 
t^z/v„  ** 


ap’ 


(3.1) 


where  f0  is  unperturbed  (equilibrium)  electron  distribution  function. 

The  Integration  must  be  done  along  the  past  history  (x'  ,p ' ,t ' )  from  the 
injection  time  t^n  "  t  -  z/v„  and  the  position  z ^  ■  0  to  the  current  time 

t  and  position  z  of  the  unperturbed  characteristic.  The  characteristic 
-  essentially  an  unperturbed  "particle  trajectory"  -  is  simply  a  helical 


trajectory  of  an  electron  in  a  uniform  magnetic  field  as  depicted  in 


Fig.  3.1  with 


-  t  +  2  ■  ~  - 


Ml 


(3.2) 


♦  '  m  $  +  u 


z  -  z 


C  V„ 

where  u>c  =  =  eHg/mc)  is  the  relativistic  cyclotron  frequency. 

The  unperturbed  electron  distribution  function  f^  is  an  arbitrary 

function  of  invariants  of  the  characteristic.  For  a  hollow  electron  beam 

they  can  be  represented  by  three  independent  invariants,  p„  pA 

» 

and  8UC**  t*'at  f0  £  where  Rg  is  the  guiding  center  radius. 


The  R  -dependency  is  essential  to  represent  an  inhomogeneous  plasma  such  as 
G 


a  finite  electron  beam.  Writing  that 


P'”  ; 


(3.3) 


one  can  cast  (3.1)  for  TE“  mode  into  the  form, 

no 


f1(x,p.t) 


I  ,  l|  y  ^Pa.  xlmcy  3p„  mcy  3p  /  w 
>t-z/v„  VL-  \  /J  <3. 


IS 


iE'  +  H’  cos©  +  H'  sinf<(i,-C^"|^--r^ 

x  mcY  z  mcy  Y  ^ /  Jnc  dRC  J 


(3.4) 
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Using  the  ansatz  given  by  (2.1),  one  can  write  (3.3)  (see  Fig.  3.1) 


f 


1 


e  f  1  .  dt* 
Jt-Z/Vu 

-  F'sinX'cose 


e“lutf  ^FjsinX’cos*’ 

+  F^cosX’isin($'  -  0)J 


(3.5) 


where 


F.  =  f*  HdS  ffo  .  -rj,  k„  ( P-i.  8f0  _  £u_  M 

1"  ykA  xk^mcy  3p„  mcy  3p_,_  y 

rl  ;  [  PI  m/ C-  _  Ml  hj,  P II  \  1  _ 0 

F2  -  ^Ey  Hx  k,  mcy  )  nc  3Rg 

-  p  t  3fn 

F*  =  H*  1  0 

3  "  z  mcY  Qc  3Rg 


and 


X’  =  X  +  a^cos^’ 

X  s  kx  (|  +  Rgcos0) 


2  +  aGcos0 


(3.6) 


(3.7) 


<aL  5  k+rL  *  aG  "  k*V 

Later,  we  will  see  that  E^-term  of  Fj^  in  (3.6)  represents  an  azimuthal 
bunching  in  the  phase  space  leading  to  a  cyclotron  maser  instability  and 
IT-term,  which  vanishes  if  fQ  has  an  isotropic  distribution  in  p„  and  px, 
represents  an  axial  bunching  (resulting  in  an  azimuthal  bunching  too)  lead¬ 
ing  to  Weibel  instability.  F^  and  F^  vanish  if  fQ  is  independent  of  Rg 
(homogeneous  plasma)  and  thus  represent  geometrical  effects  of  the  finite 
beam.  These  alone  do  not  lead  to  an  instability  but  give  substantial 
modification  to  the  behavior  of  CMI  and  Veibel,  particularly,  at  off- 
resonant  behavior. 

The  trajectory  integral  (3.5)  can  be  put  into  a  more  manageable  form 
by  harmonic  expansion  using  the  formula. 


3.3 


and,  then. 
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Then,  the  integral  (3.5)  can  be  written  as 
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(3.11) 


The  fora  of  integral  given  by  (3.9)  is  a  typical  hysteresis  integral  with 
G(z  -  z’)  being  a  kernel  ~  Green’s  function.  The  induced  current  corres- 


ponding  to  (3.9)  is  the  source  ten  in  (2.4)  making  it  an  lntegro-differen- 
tial  equation.  The  Laplace  transformation  of  (3.9)  is  given  by  convolution 


theorem  as 


where 
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F^(k),  F^k),  ^(k)  are  the  Laplace  transformations  of  F^,  F^,  F^,  respec¬ 

tively,  and,  using  (3.6)  and(3.11), 
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Notice  that  the  Interaction  is  greatly  enhanced  when  the  Laplace  transforma¬ 
tion  of  the  Green’s  function,  G(k) ,  has  poles,  i.e.,  flg(k)  *  0.  However, 
exact  locations  of  resonant  interaction  can  be  determined  from  the  full 
Maxwell-Vlasov  equation  after  we  calculated  the  Induced  current  corres¬ 
ponding  to  (3.12)  and  inserted  into  (2.8).  In  the  next  section  we  will 
calculate  the  Induced  current. 


4.1 


IV.  SQURC£  TERM-INDUCED  CURRENT 

The  induced  current  corresponding  to  the  perturbed  electron  distribu¬ 


tion  function  f ^  is  given  by 


J(x,t)  -  Ne  J*  d^p  ^  ^(x.p.t) 


(4.1) 


Where  N  is  the  number  of  electrons  per  unit  length  (the  beam  function  f^  is 
normalized  to  be  one  per  unit  length).  The  relevant  component  of  current 
is  Jy(z)  defined  by  (2.5)  and  its  Laplace  transformation  can  be  written, 
with  the  perturbed  electron  distribution  function,  (3.12)  -  (3.14),  as 
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where  G  and  S  are  defined  by  (3.8)  and 
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x  =  X  +  a^cos^ 
(X  -  SI  +  aGcos0) 


(4.3) 


First,  phase  angle  (♦)  integration  can  be  done  with  the  aid  of  formula, 

h  if"*1’*  + 


and  thus 

/  d4cos$e*8^sinx  *  -  iSJg(a^)CB  •  (4.4) 

Using  (4.4),  one  obtains  from  (4.2), 
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Next,  one  can  do  6-integration  using  the  formula 
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which  can  be  easily  proved  by  expanding  the  exponential  in  terms  of  Bessel 
functions , 
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Note  that  the  F3(k)-term  has  dropped  out  in  (4.6).  This  is  a  unique  feature 

of  a  linearly  polarized  mode  such  as  TE°  mode  where  left  and  right  heliclty 

no 

contribution  in  F3(k)-term  cancel  out. 

The  remaining  integrals  in  (4.6)  require  a  more  specific  form  of  the 
beam  function  fQ(p„,px,RG) .  However,  it  is  not  desirable  to  specialize  in 
the  beam  function  at  this  stage  and  one  useful  trick  to  avoid  this  problem 
is  to  use  the  following  Identity  relation.  For  an  arbitrary  function 

f0^P"»p**RG^  one  can  write  as 


where 
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4.3 


is  a  "cold  bean"  distribution  function  with  an  infinitely  thin  ring  guiding 
center  R^.  With  the  aid  of  this  Identity  relation,  one  can  carry  out  the 
Integration  over  Rg,  p„,  using  the  "cold  bean"  distribution  function 
and  afterward  one  can  Integrate  over  the  actual  beam  distribution  function 
f0(p«»Px»R°)  “  "statistical  average".  If  the  beam  is  actually  "cold",  the 
latter  step  ("statistical  average")  is  trivial  and  not  necessary. 

This  observation  allows  us  to  carry  out  the  integrations  in  (4.6)  for 


the  "cold  beam"  function  fg  and  postpone  the  "statistical  average"  with 
the  actual  beam  function  fg  afterward.  Understanding  this,  Eq.  (4.6)  with 
"cold  beam"  function  fg  becomes,  after  now  a  trivial  Rg- integration  with 


6-function, 
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The  remaining  momentum  Integrations  are  again  trivial  by  integration 
by  parts, 
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we 


obtain  the  final  expression  for  the  source  term  £(k)  given  by  (2.9)  as 
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If  the  beam  function  is  different  from  the  "cold  beam"  defined  by  (4.8), 
The  "statistical  average"  of  the  type  (4.7); 

<P(k)>  =  J*  2irR^dRgdp°2irp°dp°f0(p°,p°,R°)p(k)  (4.16) 

should  be  done  at  this  stage. 


5.1 


V.  GAIN  -  COLD  BEAM 

In  this  section  we  will  solve  the  Maxwell-Vlasov  equation  given  by 
(2.7),  (2.8),  and  (4.11)  for  a  cold  bean.  Using  (2.7),  one  can  write  the 


source  term  (4.14)  as 
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Substituting  (5.1)  into  (2.8),  one  obtains 
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Equation  (5.3)  gives  the  Laplace  transformed  field  E  (k)  in  terms  of  the 

—  i  d¥  ^ 

boundary  conditions  E  (0)  and  H  (0)  ■  -tt-  (0) .  The  fields  in  coordinate 

y  X  1K||  GZ 

space  can  be  easily  obtained  by  inverse  Laplace  transformation  using  the 


Bromwitz  formula 
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where  c  Is  a  sufficiently  large  positive  so  that  all  the  poles  in  F(k)  oust 


fall  above  the  integral  path.  When  F(k)  has  only  isolated  poles,  one  can 
close  the  path  by  an  upper  semi-infinite  half  circle  using  the  Jordan's 
lemma  and  the  integration  is  reduced  to  a  simple  residue  counting  at  the 
poles ; 
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In  order  to  apply  (5.6)  to  (5.3),  recall  that  Nq  is  small  (y  10  for  a  typi¬ 
cal  MIG  electron  gun)  and  only  the  term  with  ft^(k)  =  wyq  -  kp„/m  -  sflc  e  0 
contribute  substantially  in  (5.4).  In  practice,  for  a  strong  magnetic  field 


(fi  ;  large) ,  only  one  specific  harmonic  number  s  can  be  made  resonant 
c 


[JT  (k)  s  0J.  In  this  case,  (5.3)  can  be  rationalized  with  quartic  poly¬ 


nomials  as 
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where  d(k),  n^k) ,  and  n2(k)  are  obtained  from  D(k)  ,  Nj(k) ,  and  N2(k)  by 


multiplying  n^2(k).  Since  (5.7)  is  now  written  in  terms  of  rationalized 


function,  the  inverse  Laplace  transformation  (5.6)  is  trivial 
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where  k^'s  are  poles  given  by 


d(ki)  -  0 


(5.9) 


and 


d’(k)  =  ^  d(k) 


Now  we  have  completely  determined  the  field  Ey(z)  for  z  >  0  and  other 


components  of  fields  by  (2.3)  in  terms  of  the  boundary  values  E  (0)  and 


1  V 

H  (0)  •  -tt— -r^HO).  E9*  (5.8)  shows  that  the  fields  are  consisted  of 

4-propagating  modes  for  a  cold  beam  and  a  specific  harmonic  number  with 
their  relative  amplitudes  are  completely  determined.  This  is  extremely 
important  in  studying  a  backward  wave  interaction  where  all  the  4-modes 
are  important  not  like  as  the  forward  interaction  where  one  mode  eventually 
dominates  over  other  modes.  Note  that  the  condition  (5.9)  or  equivalently 
D(k)  »  0  determining  poles  is  nothing  more  than  a  usual  complex  dispersion 
relation  and  the  residues  at  the  poles  determine  the  relative  amplitudes 
among  the  modes.  If  we  had  used  Fourier  transformation  -  even  though  it  is 
ill-defined  -  we  would  have  obtained  only  the  dispersion  relation  with 
their  amplitude  information  lost. 

Typically,  the  structure  of  the  poles  are  as  shown  in  Fig.  5.1.  Other 
components  of  the  field  are  obtained  from  (5.8)  using  (2.3)  as 
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The  total  power  flow  flowing  through  the  waveguide  is  given  by,  recalling 
(2.1)  and  (2.5) 


P(t) 


f8  f 

“  I  dx 

Jo  Jo 


dy  *  §-  Fe(-  E  H*) 
7  8tr  y  x 


— * . 


2  2 
c  k„  u  /c  ab 


.BelE  (z)  •  Hx(z)} 


(5.11) 


The  quantity  which  is  more  interesting  is  gain  defined  by 


C  (z)  =  P(z) /P (0)  -  RefE  (z)  •  H*(z)]/Re[E  (0)H*(0)]  (5.12) 

r  y  x  y  x 


in  the  case  of  foward  wave  interaction  or 


FIG,  5.1  Pole  structures  in  a  backward  wave  region 
in  a  forward  wave  region. 
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C  <z)  2  P(0)/P(z)  -  R  (E  (0)H*(0)]/Re[E  (z)H*(z) ] 
o  ®  y  x  y  x 

In  case  of  backward  wave.  The  gain  In  dB  Is  defined  by 

G<dB)  =  10  log10G(z)  .  (5.14) 

A  typical  gain  as  a  function  of  interaction  length  z  or  a  function  of 
frequency  u/c  is  shown  in  Fig.  5.2. 


6.1 


VI.  THREE-WAVE  ANALYSIS:  PIERCE-TYPE  FORMULATION 

In  the  previous  section,  we  have  obtained  the  desired  solution  of  one- 

dimensional  boundary  value  problem  for  a  TE^-mode  interacting  with  a 

hollow  electron  beam.  The  solution  can  be  used  to  analyze  a  forward  wave 

amplifier  or  backward  wave  oscillator/amplifier.  The  solution  shows  that, 

at  least  for  a  "cold"  beam,  the  fields  are  given  by  four  propagating  waves 

with  their  relative  amplitudes  (including  phases)  completely  determined  by 

1  dE 

the  boundary  values ,  E  (0)  and  H  (0)  ■  -rr—  -t-^-(O)  ,  at  the  beginning  of  the 

y  x  lkfi  uz 

interaction.  Typically  two  of  them  are  real  "waveguide  modes"  propagating 
in  opposite  directions  from  each  other.  Two  others  -  "Beam  modes"  -  are 
complex  conjugate  each  other,  representing  one  growing  and  the  other  decay¬ 
ing  and  occur  near  the  point  where  the  beam  line  Intersects  with  the  wave¬ 
guide  dispersion  line.  Therefore,  three  waves  are  clustered  near  the  inter¬ 


secting  point  and  the  other  one  is  at  the  point  with  the  same  u  but  opposite 


sign  of  k„. 

As  one  can  see  in  the  following,  the  waveguide  mode  propagating 
opposite  to  the  beam  is  a  nearly  free  propagating  mode  without  Interacting 
with  the  beam  significantly.  Therefore,  one  can  eliminate  this  wave  con¬ 


centrating  of  three-wave  analysis. 

In  order  to  get  some  feelings  about  the  nature  of  the  solutions,  let 
us  write  the  denominator  of  (5.7),  whose  zeros  determine  the  possible 
modes  as  in  the  dispersion  relation  (5.9),  in  terms  of  dimensionless  para¬ 
meters.  Defining  the  "phase  shift"  h  and  the  "detuning  parameter  A  as 


h  =  k  —  k„  =  k,,£ 

°  0 
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the  'Vesonance  factor"  oP  can  be  written  as 
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Then,  factorizing  the  scale  factor  from  the  denominate  d(1c)  as 
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d(k)  =  2k,,  2*-  d(k) 
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one  obtains  the  dimensionless  expression  as 
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(6.4) 


d(0  =  c(l  +  § )fig(C)  -  g3{  [1  -  AjCd  +  f)]  +  A  Qg(0  +  A 0nJ(O)  (6.5) 


where 
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and  recall  that  is  given  by  (5.2)  and  A^(i  ■  0,1,2)  by  (4.12). 

For  a  tenuous  beam,  the  "gain  parameter"  g  is  small  in  (6.5)  and  the 
solution  of  d  -  0  is  roughly  characterized  by  the  solutions  of 

5(1  +  §>ng<5>  “  0 

which  gives  two  "waveguide  modes”  5  “  0,  5  *  -2  and  two  nearly  degenerated 

"beam  modes"  flg(C)  8  0  or  C  *  i.  The  mode  with  5  *  -2  in  (6.1)  represents 

3 

the  opposite  propagating  mode  modified  by  the  beam  only  in  the  order  of  g 
which  is  negligible.  The  other  three  modes  are  given  by  small  5  v  6  *  g 
in  (6.5)  as  a  solution  of  a  simple  cubic  equation. 


6.3 


6.4 


depending  on  the  detuning  factor).  However,  for  a  backward  wave  oscillator/ 
amplifier,  all  three  of  the  relative  amplitudes  are  Important.  Oscillation 
conditions  for  a  BWO  is  given  by  an  infinite  backward  wave  gain  defined  by 
(5.13)  and  satisfy 


(6.12) 


One  can  show  that  (6.12)  can  be  satisfied  only  by  a  series  of  specific  pairs 
of  the  frequency  and  the  detuning  factor,  (u,5).  Eq.  (6.12)  clearly 
demonstrates  that  the  backward  wave  oscillation  is  the  result  of  three 
wave  interference  and  the  correct  amplitudes  are  essential  to  obtain 
correct  oscillation  conditions. 
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Vll.  VELOCITY  SPREAD  EFFECT 

So  far  we  hava  discussed  tha  "cold  beam"  case  defined  by  (4.8).  If 
the  beam  is  general  fg(p„,  px»  Rg),  one  must  do  the  "statistical  average” 
on  the  source  term  as  (4.13).  The  moat  sensitive  velocity  effect  comas  in 
through  the  resonance  term  R?(k)  =  urrn  -  k  ^  -  sfl  and  thus  the  spread 

8  U  IQ  C 

in  the  longitudinal  momentum  is  most  Important  (spread  in  energy  Yq  is 
usually  small). 

In  this  section,  ve  will  discuss  only  this  most  Important  velocity 
spread  effect  in  the  longitudinal  momentum.  A  velocity  spread  from  a 
typical  MiG-type  electron  gun  may  be  described  by  a  generalized  Lorentian 
distribution  of  order  n 
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The  Integral  <7. 6)  can  ba  dona  most  easily  by  a  contour  Integral  over  a 
closed  contour  with  upper  or  lower  semi-infinite  circle.  Then,  the  Integral 
is  reduced  to  find  the  poles  and  their  residues  of  the  integrand.  The 
poles  are  given  by  the  roots  of 

n2°  +  1  -  0  (7.7) 

and  whatever  poles  contained  in  F(p„).  Sometimes  the  poles  of  F(p„)  occur 
either  in  the  upper  half  plane  or  lower  half  plane  depending  on  the  para¬ 
meters  but  not  in  both  planes  simultaneously.  In  this  case  one  can  close 
the  contour  to  avoid  these  poles  and  the  relevant  poles  are  entirely  from 
(7.7),  which  are  given  by 


n,  -  exp 


[-Mi- 


(i  -  l 


n) 


(7.8) 


in  the  upper  plane  and  in  the  lower  plane.  Therefore  only  either 
(1*1,  ...  n)  or  are  included  in  the  contour.  Therefore,  the  integral 
(7.6)  is  given  by 


n  n 

<F(p„)>n  *  j  ntF(p°  +  iijAp,,)/  ^ 


if  F  is  regular  in  the  upper  plane,  or 


(7.9a) 


<F(p„)>  -  ^  n*F(p°  +  n*Ap„)/  ^  ^  n* 

i-1  i*l 


(7.9b) 


if  F  is  regular  in  the  lower  plane. 

Now  one  can  do  the  statistical  average  Integral  (4.13)  with  the  source 

term  P(k)  given  by  (4.11).  The  P(k)  has  poles  in  the  complex  logltudlnal 

momentum  plane  at  Sp  =  uyn  -  k  sfi  -  0  and  thus  pJi/m  -  (uv„  -  s tl  )/k. 

s  u  m  c  u  c 

Noting  that  k  has  negative  imaginary  part,  the  poles  are  either  in  the 
upper  plane  if  u  >  s^q/yq  5  8ulc  °r  ln  the  lower  Plane  if  u  <  suc.  Usually 

the  operating  frequency  u  is  higher  than  the  cyclotron  frequency  su  for  a 

c 

forward  amplifier  (positive  doppler  shift)  and  w  is  lower  than  su  for  a  back- 

c 

ward  wave  amplifier  (negative  doppler  shift).  Therefore,  in  forward  wave 
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region  (w  >  su^) ,  one  can  close  the  contour  in  the  lower  half  plane  to  avoid 
the  poles  of  P(lt)  and,  in  the  backward  wave  region,  close  in  the  upper  half 
plane.  Then  the  statistical  average  (5.13)  picks  up  poles  in  a  forward 
wave  region  and  poles  in  a  backward  wave  region.  Using  (7.9a)  or  (7.9b), 


one  obtains 


*\»  TT — ■>  A'X'  n  e  * 

<p(k)>  -  niP(k,p„  +  njAp,,)/  y.  n± 


(7.10a) 


in  a  forward  wave  region  (u  >  su  )  or 

c 


II  II 

<P(k)>  -  n1P(k,p?,  +  r^Ap,,)/  nA 


(7.10b) 


in  a  backward  wave  region  (w  >  sw^) , 


With  this  new  source  term,  (7.10a)  or  (7.10b),  Eq.  (5.4)  is  replaced 
by 

D00  -  (k2  -  k2>  -  »„  jr  [(4  -  k2W  Aj  <-lj  >+si)cA1<  ±5  >  +  A01 

S*»- « L '  '  S  8  -I 


Nx(k)  *■  k  + 


N2(k)  -  k„ 


CD 

„  VTk“-A  <-i-. +<!&•>  A1 

0  2J  m2  2  n°2  00  M 

S«-w  s  s  -J 


(7.11) 


in  a  backward  wave  region,  or  the  same  expression  as  (7.12)  but  is 
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replaced  by  In  a  forward  wave  region. 

A  major  qualitative  difference  of  (7.11)  from  (5.4)  Is  that,  when  they 
are  rationalized  as  (5.8),  d(k)  is  no  more  a  quartic  polynomial  If  n  £  2. 
This  means  that  the  number  of  waves  are  more  than  four.  For  n-1 ,  d(k)  is 
still  a  quartic  polynomial  and  described  by  four-waves  as  In  the  cold  beam 
case.  However,  as  one  can  see  that  an  ordinary  Lorentian  distribution 
(n*l)  has  too  long  a  tail  and  may  be  unrealistic.  The  case  with  n**2  is 
seen  to  be  much  more  realistic  and  then  the  number  of  waves  are  six.  In 
general,  the  possible  number  of  modes  are  2  +  2n.  Here  we  saw  that  the 
number  of  waves  required  to  account  for  a  velocity  spread  effect  is  highly 
dependent  on  the  shape  of  velocity  distribution  function. 


VIII.  DISCUSSION 


In  this  report  we  have  presented  a  linear  theory  to  analyze  gyrotron 
BWO.  This  has  been  possible  because  our  formulation  allows  that  the  axial 
behavior  of  the  rf  field  and  the  electron  beam  is  completely  determined 
self-consistent ly  by  the  boundary  values.  This  type  of  solving  the  problem 
as  a  boundary  value  problem  is  particularly  crucial  for  analyzing  BWO 
because,  in  BWO,  the  electron  beam  Itself  plays  the  role  of  self-feedback. 

It  has  turned  out  to  be  that  the  infinite  backward  gain  condition  in 
a  BWO  can  be  met  by  an  interference  among  three  backward  propagating  modes 
with  appropriate  amplitudes  (including  phases)  and  growing  and  decaying 
factors.  If  the  conditions  are  right  with  the  detuning  factor  and  inter¬ 
action  length,  rf  fields  can  grow  out  of  a  noise.  A  typical  backward  wave 
oscillation  condition  as  a  function  of  the  detuning  frequencies  and  inter¬ 
action  lengths  are  shown  in  Fig.  8.1.  There  are  a  series  of  oscillation 
points  with  different  detuning  factors  and  interacting  lengths.  Typically 
the  smaller  the  detuning  factor,  the  longer  the  interacting  length.  On  the 
actual  device,  it  is  interesting  to  see  whether  one  can  see  all  these  modes 
or  only  the  first  mode  (with  the  shortest  interaction  length)  dominates  and 
blocks  out  the  remaining  modes.  The  experimental  design  has  a  provision  to 
test  this.  With  a  reasonable  interaction  length  less  than  15  cm,  tuning 
over  almost  an  octave  in  Ku-band  seems  to  be  possible.  The  velocity  spread 
effect  does  not  deteriorate  much  up  to  10Z  or  so,  worse  than  the  gyrotron 
oscillator  but  better  than  the  forward  wave  amplifier.  Another  interesting 
effect  of  velocity  spread  is  that,  when  it  is  above  a  certain  critical 
value,  it  stabilizes  the  oscillation  -  acts  as  a  backward  wave  amplifier 
rather  than  oscillator.  This  opens  up  a  possibility  to  operate  a  gyrotron 
BWO  in  a  backward  wave  amplifier.  Certainly  the  required  interaction  length 


FIG.  8.1  Oscillation  conditions  as  a  function  of  detuning 
frequenc;  and  interaction  length  for  different 
current  of  the  beam. 
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is  large  (typically  30-40  cm)  and  the  efficiency  may  be  low  due  to  the 
large  velocity  spread. 

In  conclusion,  the  analysis  shows  that  a  several  tens  killowatt  gyro- 
tron  BWO  may  be  possible  in  a  Ku-band  tunable  up  to  almost  an  octave  by 
either  magnetically,  electrically,  or  both.  The  tunability  may  find  an 
ideal  application  in  studying  the  rf  heating  in  fusion  devices. 
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